Rows of optical vortices from elliptically perturbing a high-order beam 
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An optical vortex (phase singularity) with a high topological strength resides on the axis of a high-order 
light beam. The breakup of this vortex under elliptic perturbation into a straight row of unit strength 
vortices is described. This behavior is studied in helical Ince-Gauss beams and astigmatic, generalized 
Hermite-Laguerre-Gauss beams, which are perturbations of Laguerre-Gauss beams. Approximations of these 
beams are derived for small perturbation, in which a neighborhood of the axis can be approximated by a 
polynomial in the complex plane: a Chebyshev polynomial for Ince-Gauss beams, and a Hermite polynomial 
for astigmatic beams. © 2008 Optical Society of America 

OCIS codes: 140.3300 (laser beam shaping), 260.2110 (electromagnetic theory), 999.9999 (optical vortices) 



The light beams most often studied have cylindrical sym- 
metry - their intensity is invariant with rotation about 
the optic axis. When such beams carry orbital angular 
momentum,^ such as Laguerre-Gauss (LG)^ and Bessel 
beams, they possess an axial optical vortex (phase sin- 
gularity),''^ where the intensity vanishes, and about which 
the phase changes by 27r^. This integer vortex strength 
£ equals the beam's orbital angular momentum. For a 
high-order beam > 1), the axial singularity is un- 
stable to perturbation. A simple example is the addition 
of a small-amplitude cylindrical beam without a vortex 

= 0), whence the axial vortex unfolds into \^\ unit 
strength vortices which, in a transverse plane, are equally 
spaced on a circle centered on the axis. Thus rotational 
symmetry is broken from continuous to discrete. 

Here, I want to draw attention to a natural class of 
perturbations, where the high-order vortex breaks up 
into a straight row of same sign unit strength vortices, 
destroying rotational symmetry. At least three exam- 
ples exhibit this: helical Ince-Gauss (IG) beams, ^'^ re- 
lated to gaussian modes separated in elliptic coordinates; 
their propagation-invariant analogue, Mathieu beams;^ 
and Hermite-Laguerre-Gauss (HLG) beams,* generated 
by astigmatic transformations of LG and Hermite-Gauss 
(HG) beams. IG and HLG beams are perturbations of 
LG beams; Mathieu beams, of Bessel beams. 

The unfolding of high-order vortices into rows under 
certain transformations is important in the understand- 
ing of the physical structure of vortex cores. In par- 
ticular, in any synthesis of LG beams from HG modes 
with mode converters. any misalignment will lead to a 
row of closely spaced vortices, as will be described. 

Two assumptions are made without loss of generality: 
the description is confined to the waist plane, and i is 
positive. A normalized LG beam, of mode order = 
£ -I- 2p, is therefore represented^ 



[2''+V!]^/^i?''exp(i£0) 



2R^ 



(1) 



where Li^ is an associated Laguerre polynomial, R, 



are polar coordinates, and Wq is the waist width. It is 
proportional to exp{i£(j>) = {x + iyY near the axis, 
like any complex high-order beam. 

A beam's nodal structure may be understood in terms 
of crossings of the zero contours of its real and imaginary 
parts. For R^ exp{i£(j)) , the real and imaginary contour 
pattern consists of 2£ alternating, regularly spaced ra- 
dial lines, reflecting the £-fold rotational symmetry of 
the phase singularity (see Fig. ^). Adding a small real 
constant a > (approximating a vortcx-free cylindrical 
beam) only affects the real contours (Fig. QJd) , leading to 
vortices at positions — a^/^ exp(27ri7i/£) for n = 1, 

This contrasts with elliptically perturbed beams, 
where the high-order vortex breaks into a row of £ 
vortices of strength -1-1. These beams are conjugation- 
symmetric about the x-axis (real (imaginary) part 
(anti)symmetric), with an imaginary zero at y = 0; vor- 
tices occur where the real zero contours cross this line. 
These vortices have the same sign, which appears to vi- 
olate the sign rule,^^ since they all lie on the same imag- 
inary zero contour. However, between each pair of real 
contours, another imaginary contour crosses the x-axis - 
there is a phase saddle between adjacent vortices, which 
must therefore have the same sign by the extended sign 
rule.^^ Examples are shown in Fig.^ and d. 

The presence of these saddles illustrates a general fea- 
ture: upon perturbation, a high-order vortex unfolds, not 
only to £ unit strength vortices, but also £ — I saddle 
points. This is because the Poincare index of phase gra- 
dient (current) must be conserved; a vortex, of whatever 
strength, is an index circulation, so must be balanced 
by £ — 1 index —1 phase saddle points. In perturba- 
tion by a constant, these saddles remain degenerate at 
the origin. These observations will now be justified for 
helical IG beams, Mathieu beams, and HLG beams. 

For the discussion of IG beams, the notation of Ref. 
5 will mostly be adopted. In the waist plane, the real 
and imaginary parts of helical IG beams are standing 
IG modes (real part symmetric, imaginary part antisym- 
metric). Each mode is a product of Ince polynomials^' 
in elliptic coordinates u, v; contours of constant u (v) are 
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Fig. 1. Intensity, real (solid) and imaginary (dashed) 
zero contours for beams in the waist plane: (a) LG beam 
V'LG,4,o; (b) cylindrically perturbed ^/'lg,4,o + O.2'0lg,o,o; 

(c) helical IG beam i/^ig 4 4' (d) HLG beam i/'hlg 4 o- (^) 
is the unperturbed beam corresponding to the others. In 
the elliptically perturbed cases (c) and (d), the order 4 
vortex has unfolded into a row of four strength 1 vortices, 
interspersed by saddles with the same phase i. 



confocal ellipses (hyperbolas) with foci at y = 0,x ^ ±/o 
(u + iv = arccosh((x + iy)//o)). The dimensionless ellip- 
ticity parameter e = 2/q/wq is used. When e = /o = 0, 
IG beams arc LG beams (elliptic coordinates become po- 
lar), and as /o, e — > cx) they become cartesian HG beams. 
A helical IG beam may therefore be treated as an ellip- 
tic perturbation of a LG beam for /o,e > 0. Real IG 
modes are written as a product of exp{—R^/wQ) and a 
pair of Ince polynomials in v and m; they therefore have 
zeros on ellipses and hyperbolas. The helical IG beam of 
ellipticity e corresponding to ip-LG,i,p is therefore^ 



V'lG.w/ = [^c^C^(iM,e)C^(w,e) 

+ ^^^5^(m, e)5^(v, e)] exp(-i?V«;g), (2) 

where the A denote appropriate normalization constants 
(found in Ref. 17 Eq. (3.14)). The Ince polynomials 
Cf^,Sff are finite trigonometric series solving Incc's 
equation;^^ Cf^ is a sum of cosines, 5^ of sines. In par- 
ticular, C^{iu,e) is a sum of cosh terms, S^{iu,e) of 
i sinh terms, so the antisymmetric part is imaginary. 

When e > 0, the pattern of real and imaginary zeros 
of a high-order vortex therefore unfolds into a system 
of confocal hyperbolae, as shown in Fig. given by 
the zeros of Cf^{v, e) and Sj^{v, e). These zeros {£ for C, 
£ — 1 for S) alternate, since both functions are eigen- 



functions of Ince's equation, which is of Sturm-Liouville 
type: between adjacent zeros of S^{v,e), Cf^{v,e) must 
have a zero."'^* Thus helical IG beams possess a straight 
row of alternating equal-sign zeros and saddles between 
their foci, with positions corresponding to the zeros of 
C^(arccos(x7/o), e) and 5^(arccos(x//o), e). 

The transition from LG to helical IG beams has been 
described as a perturbation; to sec how the high-order 
vortex becomes a row, it is instructive to study the limit 
of small e. As e ^ 0, Ince's equation tends to the 
equation of simple harmonic motion, and C^(w,e) — > 
cos{£v), Sj^{v, e) sin{£v). Thus, for small e = 2/Q/wg, 

^IG Af £ ~ ^ (cos(fo) cosh(£u) + i sin{£v) sinh(£u)) 

= ^cosh(£(u + iv)) = ATtiix + iy)/fo)), (3) 

where Tg is a Chebyshev polynomial of the first kind^^ 
(in the limit, the normalizations of C and S are equal^^). 
A neighborhood of the origin, of order /g, may there- 
fore be approximated by a polynomial in the complex 
plane, whose zeros, which scale with /o, are all on the 
real axis. A consequence of this representation is that all 
zeros have sign -1-1 (they are not poles), and their cores 
are isotropic. This isotropy is somewhat surprising, as 
one might expect the singularity phase structure to be 
squeezed as the row contracts; the isotropy is ensured by 
the complex analytic approximation.^" 

The argument for Mathieu beams is almost identical. 
The vortex (saddle) positions on the row is given by the 
zeros of (anti)symmetric Mathieu functions, which have 
the same limiting behavior as Ince polynomials. 

HLG beams are the final example. They occur in ex- 
periments when a HG or LG beam undergoes an astig- 
matic transformation,® due, for instance, to a cylindrical 
lens or variable-phase mode converter, the perturba- 
tion parameter being the orientation angle a of the lens. 
(In the language of Ref. 11, 6 = tt/2 and (f> = a.) They 
can be mathematically understood using the analogy be- 
tween gaussian beams in the waist plane and the quan- 
tum 2D harmonic oscillator:^^ HG states correspond to 
linear, LG to circular, and HLG to elliptic orbits. Using 
Schwinger's analogy between the 2D harmonic oscilla- 
tor and quantum spin,^^ HLG beams may be written as 
sums of LG beams whose coefficients are spin rotation 
matrix elements (Wigner d-functions). Therefore, the LG 
beam iphCi.p is perturbed by a to 



N/2 

HLG/,p = ^ ' 

m=-N/2 



,N/2 

'e/2,m 



(2a)(-l)"-|"IVLG,2,n,Ar/2-|rr«| 

(4) 

where N ~ £ + 2p and the d-function d^^, is given 
in Ref. 23 Eq. (3.8.33). Note that for a = 'n/A, -0hlg 
is a HG beam, and an alternative representation of 
Eq. gl) has a replaced by 7r/4 - a, '0LG,2m,Ar/2-|m| by 
V'HG,JV/2+m,Af/2-m ^ud different phase factors. This sec- 
ond form is equivalent to that given in Ref. 8 Eq. (8) in 
terms of Jacobi polynomials, and for a = 7r/4, in Ref. 21 
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Eq. (19) using d-functions. Rows of vortices have been 
experimentally observed for these beams.*' 

Unlike IG beams, HLG beams cannot be expressed in 
terms of a separable coordinate system (since only carte- 
sian (HG), polar (LG) and elliptic (IG) are possible^^). 
Thus, vortex breakup into a row cannot be justified using 
Sturm-Liouvillc analysis. Furthermore, for large enough 
a, the negative vortices can appear on the row. which 
ultimately vanishes (when a = 7r/4, the beam is HG). 
However, the high-order vortex of iphG does break into a 
row of vortices and saddles for small a, and in fact, may 
be approximated by a complex analytic function. 

Approximation of Eq. Q requires that i?, a sa 0. To 
first term in i?, Eq. is 

^ [{\e\+p')\]'/^x + ^sign{e')y)\''\ 



(7rp'!)i/2|f 



(£ may be positive or negative, since all signs appear in 
Eq. I0J). For small a, c?^, m(2'^) can be approximated 



C™(2«) 



U + m'y.{j - m)\ 



(j — m'y.{j + m)! 



1/2 



(-ir 



(m' 



(6) 

when m' > 0,m' > m (when m > m', then m^m! are 
reversed and there is no sign factor). Substituting these 
into Eq. (0J, and using the leading-order scaling / a ^ 
constant, it can be shown that 



HLG,£,p 



T:p\ 



ly 



(7) 



where He is the Hcrmitc polynomial. That is. for 
small a and R, a HLG beam is proportional to a Hermite 
polynomial in the variable (x + iy)/2wQ^/a; this is again 
a complex analytic approximation, so the £ zeros have 
strength one and are isotropic. The vortices lie in a row 
since Hermite polynomials have real zeros. 

I have described the phenomenon of the breakup of a 
high-order phase singularity in a rotationally symmet- 
ric beam into a straight row of equal sign, unit strength 
singularities, and gave explicit examples for helical Ince- 
Gauss and generalized Hermite-Lagucrre-Gauss (astig- 
matic) beams. In each case, the unfolded row of £ vortices 
scales by the square root of the dimensionless perturba- 
tion parameter (e or a), as this tends to zero. In this 
limit, the transverse neighborhood of the beam axis is 
approximated by a polynomial depending on the com- 
plex variable x + iy, with real zeros, automatically im- 
plying a straight row of vortices with equal sign. The 
intensity near the row is also proportional to the pertur- 
bation strength (given explicitly in Eq. (O). 

Helical IG and HLG beams are distinct, since Eqs. JSJ 
and involve different polynomials; each row's rela- 
tive vortex spacings are different. It is likely that elliptic 
perturbations can be characterized using group theory,^'' 
but this is outside the scope of the present work. 

In a paraxial beam, an optical vortex is a trans- 
verse solution of Laplace's equation, so it has the form 



{x + iyY + a{x — iyY , with \a\ < 1.^^ Cylindrical symme- 
try implies an axial vortex is isotropic (a = 0). Isotropy 
of the unfolded vortices implies that the perturbation, 
to leading order, only involves cylindrical beams of the 
same sign. An anisotropic, high-order vortex thus will 
unfold into anisotropic vortices, and not approximable 
by a function of a; -f- iy. Of course, there are many other 
ways an order £ LG beam may be perturbed, involving 
beams of the same mode order. However, it is quite dif- 
ficult to construct superpositions in this way leading to 
vortex rows; nevertheless, the elliptic perturbations de- 
scribed here are physically natural. 
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